Inelastic Collisions of a Fermi Gas in the BEC-BCS Crossover 
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We report the measurement of inelastic three-body and two-body collisional decay rates for a 
two-component Fermi gas of ®Li, which are highly suppressed by the Pauli exclusion principle. Our 
measurements are made in the BEC-BCS crossover regime, near the two-body coUisional (Feshbach) 
resonance. At high temperature (energy) the data shows a dominant three-body decay process, 
which is studied as a function of bias magnetic field. At low energy, the data shows a coexistence 
of two-body and three-body decay processes near and below the Feshbach resonance. Below reso- 
nance, the observed two-body inelastic decay can arise from molecule-atom and molecule-molecule 
collisions. We suggest that at and above resonance, an effective two-body decay rate arises from 
collisions between atoms and correlated (Cooper) pairs that can exist at sufficiently low temperature. 



Quantum statistics dramatically affects the inelastic 
collision rates that determine the lifetime of cold atomic 
gases. In an inelastic three-body collision, two of the 
colliding atoms decay to a bound molecular state, releas- 
ing energy. Interactions between atoms can be strongly 
enhanced by tuning a bias magnetic field near a colli- 
sional (Feshbach) resonance In a Bose gas, this 
enhancement is accompanied by an inelastic collision rate 
that increases by two or three orders of magnitude com- 
pared to that obtained away from resonance [s'l, and a 
correspondingly short lifetime of just a few ms at typ- 
ical atomic densities. In contrast, for a Fermi gas in a 
mixture of one or two different spin states, the proba- 
bility of three atoms colliding is highly suppressed by 
the Pauli exclusion principle. The lifetime of the cloud 
is on the order of 0.1 s for fermionic [1, and 
50 s for ^Li [^, The long lifetime of Fermi gases 
is essential to the study of strongly interacting Fermi 
gases 0, P , which offers unprecedented opportunities to 
test nonperturbative theoretical techniques that apply to 
exotic systems ranging from high temperature supercon- 
ductors to nuclear matter. Determination of the inelas- 
tic collision rate coefficients in the strongly interacting 
regime of a Fermi gas provides new tests of few-body 

theories [13, El, [H, El, H, El, H [13, El ■ 

In this Letter we report on the precision measurement 
of three-body inelastic collision rate constants for 
an ultracold two-component Fermi gas in the BEC-BCS 
crossover regime near a Feshbach resonance. We also ob- 
serve two-body inelastic decay below the Feshbach res- 
onance, which arises from molecules [3, 16 1. From the 
data, we estimate the corresponding rate constants K2. 
Finally, we observe two-body decay at and just above the 
Feshbach resonance. We suggest that this process arises 
from correlated pairs, which is a many-body effect. We 
load a Fermi gas from a single beam CO2 laser trap into a 
CO2 laser standing wave that is formed by the incoming 
and retro-reflected beam. The standing wave produces a 
potential with a period of 5.3 /im that is four times deeper 
than that of the single beam trap and tightly confining in 
the axial direction (along the standing wave). The cor- 



responding atomic density is up to lO^^'/cm'^, ~ 20 times 
higher than that obtained in the single optical trap. This 
dramatically increases the inelastic collision rates, mak- 
ing precise measurement of the rate constants feasible. 

For two-component Fermi gases, three-body inelastic 
collisions arise in the BEC-BCS crossover for processes 
of the form F + F + F' ~* F+{FF'), where F and F' are 
fermions in different states and {FF') is a bound molec- 
ular state. On the BEC side of the Feshbach resonance 
where the scattering length a > 0, the three-body decay 
rate is predicted to scale as 
side (a < 0), it should scale as 
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15|, while on the BCS 
|a|2 ''55 [2]. By contrast, 
two-body inelastic collisions can arise from the decay of 
real molecules, which exist on the BEC side. These pro- 
cesses take the form either {FF') + F ^ {FF')^ + F or 
{FF') + {FF') {FF')^ + {FF'), where (FF')_ is a 
deeply bound molecular state. The theory predicts that 
the decay rate scales as a~^'^^ for atom-molecule colli- 
sions or a~^'^^ for molecule-molecule collisions [ll|. 

In the experiments, a sample of ^Li atoms in a 50-50 
mixture of the two lowest hyperfine states is loaded into 
a CO2 laser trap with a bias magnetic field of 840 G, 
where the two states are strongly interacting. Evapora- 
tive cooling is performed to lower the temperature of the 
sample [1|. The magnetic field is then changed in 0.8 sec- 
onds to a final magnetic field where we perform the mea- 
surement. Subsequently, the gas is adiabatically loaded 
into a CO2 laser standing wave by slowly turning on the 
retro-reflected CO2 laser beam. A quasi-two-dimensional 
Fermi gas is then formed and absorption images are taken 
at various times after the formation of the 2-D system to 
determine the inelastic decay rate. 

At the final optical trap depth, the measured trap os- 
cillation frequencies in the standing wave are uj± = 27r x 
3250 Hz in the transverse directions and = 27r x 83.5 
kHz in the axial direction. The corresponding frequen- 
cies in the single beam trap are uj± = 2tt x 1650 Hz and 
a;^ = 27r X 56 Hz, respectively. Our measurements in- 
dicate very good standing wave alignment, as the trans- 
verse frequency is nearly twice that of the single beam 
trap, as expected. 
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FIG. 1: Atom number versus time. Data were taken at 834 
G and E±/Ep± = 1.8. A'^ is total atom number and No is 
initial atom number in the observed region of the cloud. Blue 
dots: Experimental data; Red solid curve: Three-body decay 
fit; Green dashed line: Two-body decay fit. 



FIG. 2: Three-body inelastic collision rate coefficient K3 ver- 
sus atomic density for iJx/i^Fx = 1.8. Blue dots: Experimen- 
tal data. Error bars indicate statistical errors; Red dashed 
line: Fit to the data whh K3 = (8.44 ± 1.04) x 10~^*^cm*'/s. 



The total energy of the gas obeys the virial theo- 
rem [19I when the bias magnetic field is tuned to a 
broad Feshbach resonance, where the Fermi gas is uni- 
tary. Since the trap depth is large compared to the energy 
of the cloud, the confining potential U is approximately 
harmonic. Then the total energy is E = 2{U) — Ez+E±, 
where E^ is the axial energy and E± is the transverse en- 
ergy, referred to the trap minimum. We determine only 
the transverse energy E± — 2muj'^{x'^) , by measuring the 
mean square transverse cloud size (x^) . For reference, the 
transverse energy for the ground state of an ideal two di- 
mensional Fermi gas is Ej± = ^Ep±, where Ep± is the 

1/2 

transverse Fermi energy, Ep± — huj±Ns . Here m is 
atomic mass of ^Li and Ng is the total atom number in 
one site. For our experiments in the unitary gas, we mea- 
sure E±/Ep±_ - 1.8 with Ns = 2, 600 and Ejl/Ep± ~ 0.7 
with Ng — 1,600. If the 2D unitary gas has the same 
effective mass as the 3D case, the 2D ground state trans- 
verse energ y w ould be 2Ep±y/T+^/3 ~ 0.42 Ep±, using 
/3 = -0.60 [20|. In this case, our lowest energy would be 
significantly above the ground state value. 

In general, for magnetic fields away from resonance 
where the scattering length is finite, the total energy is 
dependent on the scattering length In this case, 

we measure the number-independent mean square trans- 
verse cloud size {x'^)/xpj_, where Xpj_ is defined by 
2mui']_x\^ = Ep±. For an ideal gas in the ground state. 
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we note that (xq) 

We measure inelastic collision rates by measuring the 
time dependence of the atom number and the radial cloud 
size. The atom number TV as a function of time is Q 



dN 
~dt 



= -TN - K2n^ (Tx ~ K^.n'' (Fx 



(1) 



where n is the atomic density. On the right side, the first 
term arises from background collisions with a density- 
independent rate F (l/F = 64 s for our trap). The second 



term arises from loss due to two-body inelastic collisions 
with a rate coefficient K2, while the third term arises 
from loss due to three-body collisions with a rate coeffi- 
cient K^. 

For the conditions of our experiments, where 
Ep^/fiLUz — 1.5, the ground axial state contains 90% 
of the atoms for an ideal Fermi ga at zero temperature. 
For simplicity, assume that the 2-D Fermi gas is primar- 
ily in the ground axial state of a single site. Then, the 
atomic density is 



n{p,z) 



2 N{z) 



1 



exp 



(2) 



for < p < (7±. Here, N{z) is atom number in the 
site at position z. a ^ is transverse width for a fit of a 
Thomas-Fermi distribution to the atomic density profile 
in the transverse directions, = (— ^)i/^ is axial width 
for the ground state (along the standing wave), and 
is the corresponding axial trap frequency. 

In our experiments, N(z) varies as a gaussian distribu- 
tion function of z with width over the whole cloud in 
the axial direction. Strictly speaking, cr_L , and lOz also 
vary with z since the depth [/(z) of the potential for a 
site at z is a Lorentzian function of z. However, we mea- 
sure a restricted part of the cloud from z — — 0.83Lz to 
z = 0.83 Lz over which lJ{z) varies less than 10%. Hence, 
to good approximation, crj^, Oz and are spatially con- 
stant. 

Integrating the atomic density over each well and then 
over the restricted region of the cloud, we obtain from 
Eq.[T] 



dt 



a2K2 2 



(3) 



where is total number of atoms in the restricted re- 
gion. Here a2 — '^^n^'^^'^ and — -^n^^. Note that 
a± {t) is a function of time since heating leads to an in- 
crease in temperature and hence the width of the cloud 
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during the atom loss process. Typically CT']_{t) and aj_{t) 
can be fit well to exponential curves, oc exp(7t). Note 
that at the highest energies used in our experiments, a 
significant fraction of atoms can occupy the first axial ex- 
cited state. If we assume a 50% fraction, the coefficient 
as is decreased by a factor 0.78, while a2 is decreased by 
a factor 0.88. These systematic corrections are smaller 
than the statistical uncertainty in our data, so we neglect 
them in our initial analysis. We then can assume that the 
axial width is time independent. 

In the first set of experiments, we have measured atom 
number as a function of time in the unitary regime at 
the Feshbach resonance (834 G), as shown in Fig.[T] The 
trap depth is set at 20% of the maximum attainable by 
reducing the laser intensity. The measured transverse 
energy of the cloud is E±/Ef± = 1.8. We observe a sig- 
nificant (> 60%) loss of the atoms in ^ 20 sec. The data 
is fit with Eq. [3l We find that a three-body decay curve 
fits the data very well while a two-body decay curve does 
not. This indicates that three-body inelastic collisions 
play a dominant role in the atom loss. 

Fig. shows the inelastic decay rate coefficient K3 as 
a function of atomic density, at the Feshbach resonance, 
for Ej_/Ep± = 1.8. The atomic density is varied by 
varying the final trap depth. Data are fit to three-body 
decay curves, from which we determine K3. A constant 
value of K3 over factor of 10 in atomic density indicates 
the atom loss is indeed a three-body decay process. By 
fitting all of the data with the same K3, we obtain = 
(8.44 ± 1.04) X 10-28cm6/s. 

We have also measured K3 as a function of magnetic 
field for {x^)/x'^j_ = 1.8, which corresponds to the trans- 
verse energy E±/Ef± = 1.8 at resonance. The fitted K3 
is plotted as a function of interaction strength l/kp±a, 
Fig. H Here kF± = {2mEFi_f''^/Ti is the two dimen- 
sional Fermi wave vector for an ideal gas at the trap cen- 
ter and a is the s-wave scattering length. By tuning the 
magnetic field from 790 G to 1200 G, we vary 1/fcp 1 a 
from 0.20 to -0.56, using the known values of a[B) [22| . 
A factor of ~ 40 decrease in is observed as the bias 
magnetic field is tuned from the BEG regime to the BGS 
regime. We fit our data on the BGS side of the Fesh- 
bach resonance with the function of = C|a|" and find 
n = 0.79 ± 0.14. The result is in significant disagree- 
ment with the theoretical prediction n — 2 Abb [l^. On 
the BEG side, increases as the magnetic field is tuned 
away from the Feshbach resonance, instead of peaking on 
the resonance. This is consistent with the experiments 
by other groups [1, H, 0] ■ 

We have repeated the measurement of atom number 
versus time, at resonance in the unitary regime, but at a 
lower energy E^/Efi. = 0.7, Fig. H Neither two-body 
decay alone nor three-body decay alone fits the data. In- 
stead, the combination of two-body and three-body decay 
fits the data well, which indicates two-body and three- 
body decays both contribute to the atom loss. 
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FIG. 3: Ks versus interaction strength l/kp±a at {x^)/x'^^ = 
1.8. Bars denote statistical error. Varying the magnetic field 
from 790 G to 1200 G changes l/kF±a from 0.20 to -0.56. 
We observe a factor of 40 change in K3, from (17.3 ± 3.2) x 
10"^*cm®/s at 790 G to (0.44 ± 0.22) x 10"^*cm®/s at 1200 
G. 

We suggest that the two-body process is related to 
correlated pairs that can exist at low energy (temper- 
ature). At higher energy, only single atoms exist while 
pairs are broken. In that case, the Fermi gas can only de- 
cay through three-body inelastic collisions of free atoms. 
By contrast, at low energy, pair-atom or pair-pair in- 
elastic collisions are possible. Therefore, both two-body 
decay and three-body decay processes can play a role in 
the atom loss. 

By measuring atom loss as a function of time at 
E^/Ef± = 0.7, we find K3 = (3.30± 1.81) x lO^^Scm^/s 
and K2 = (0.42 ± 0.16) x lO'^cm^/s. It appears that 
K3 is approximately three times smaller than that at 
E±/Ef± = 1.8. This suppression cannot arise from Pauli 
blocking, as the energetic final states are unoccupied. 

The observed scaling of K3 with transverse energy is 
consistent with the prediction of Ref. [l3|, where oc 
E for the lowest order process. We observe Ks{Efi_ = 
1.8) /K3{Ef± = 0.7) = 8.44/3.30 = 2.56, in very good 
agreement with the predicted ratio, 1.8/0.7 = 2.57. 

Although the data indicates a linear scaling of with 
energy, a decrease in can also arise from a reduc- 
tion in the number of available single atoms, due to pair 
formation. Defining / as the fraction of atoms which 
are paired, the three-body decay rate is proportional to 
(1 — f)^N^. For pair-atom collisions, a two-body rate 
would scale as /(I — f)N'^, while for pair-pair collisions, 
the corresponding rate would be proportional to pN-^. 

Using these assumptions, we can rewrite the rate con- 
stants that appear in Eq. [3] as 

K2 EE fK°,^fKl^^ + f{l-f)Kl^,. (4) 

Here K2 pa is the pair-atom inelastic collision rate coef- 
ficient and K2 pp the pair-pair inelastic collision rate 
coefficient. At E±/Ef± = 1.8, we observe pure three- 
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FIG. 4: Atom number versus time. Data were taken at 
E±/Efi_ = 0.7 in the unitary regime. Blue dots: Experimen- 
tal data; Red solid curve: Combination fit including two-body 
and three-body decay; Violet dotted line: Three-body decay 
fit; Green dashed line: Two-body decay fit. 



over a wide range of density, we are able to distinguish be- 
tween two-body and three-body processes. This method 
may provide a probe to determine the fraction of pairs or 
molecules in the Fermi gas, once the energy scaling of 
is fully established. In the unitary regime, investigation 
of the energy (or temperature [20|) dependence of K3, as 
well as the pair fraction, will be an important topic of 
future work. 
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Foundation, and the Chemical Sciences, Geosciences and 
Biosciences Division of the Office of Basic Energy Sci- 
ences, Office of Science, U.S. Department of Energy. We 
are indebted to Le Luo and Bason Clancy for help in the 
initial stages of this work. 



body decay so that / = 0. Hence we have — if 3 — 
(8.44 ± 1.04) X 10-28 cm^/s. 

If we make the extreme assumption that is indepen- 
dent of energy, then we can reinterpret the fitted values 
of K3 and K2 for E±/Ef± = 0.7 using Eq.[4]for the rate 
constants. K3 = {1 - ffK^ yields / = (30 ± 15)% and 
K2 = fK^ then requires = (1.72±1.04)x IQ-i-^cm^/s. 
As the fraction of pairs appears large, it is more likely 
that the reduction in arises at least in par t from en- 
ergy scaling, which agrees with predictions [l3|, and that 
the true fraction of pairs is smaller. 

At a magnetic field of 790 G, we first analyze the data 
to determine K3 and K2 of Eq.[3l For {x'^)/xj^j_ = 1.8, we 
find K3 = (17.3±3.2) X lO-^Scm^/s. At (x^)/xl^^ = 0.9, 
we obtain, A'3 ^ (9.35±3.06) x lO^^^cm^/s, which is con- 
sistent with the predicted linear scaling with energy [l3| . 
The corresponding two-body decay rate constants are 
if2 = at (a;2)/x|._L = 1.8 and K2 = (0.57 ± 0.22) x 
10-"cmVs at (x2)/x|,j^ = 0.9. 

If we again assume instead that of Eq.[4]is indepen- 
dent of energy, we have K° = (17.3 ± 3.2) x lO'^^cm^/s. 
Using K3 = (9.35 ± 3.06) x lO^^s^j^e/g ^j. {x^)/xl^ = 
0.9, we require the molecular fraction to be / = (19±7)%. 
Then, we obtain = (3.22 ±0.60) x lO-i-^cm^/s. Note 
that, on the BEC side, two-body inelastic collisions are 
expected to be molecule-atom or molecule-molecule, as 
predicted fl6| . The increased two-body rate arising from 
molecules on the BEC side supports our assumption that 
the two-body rate at and just above resonance arises from 
correlated pairs. In this case, a many-body theory of in- 
elastic collisions will be needed to replace the few-body 
theory that is valid far from resonance. 

Above the Feshbach resonance, we do not observe a 
two-body decay process for l/{kp±a) < —0.09, i.e., B > 
860 G. This suggests that no pairs are formed for B > 860 
G at the lowest energy E±/Ep± — 0.7 we achieve. 

By comparing the data at high energy and low energy 
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